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Cryogenic Propellant and Pressurization Modeling 


♦ Presentation Outline 

• Computational Propellant and Pressurization Program - One Dimensional 
(CPPPO) Model Development 

- Define Control Volumes 

- Various Model Information 

- Conservation Equations 

- Heat and Mass Transfer 

• CPPPO Model Validation 

- Conservation of mass and energy internal to model 

• Validated during all scenarios for heat and mass transfer 

- Comparison to existing analytic model 

• ROCETS used for design of Ares I Upper Stage 

- Flight Experiment 

• AS-203 (S-IVB stage orbital test) 

- Ground Experiment 

• NASA Glenn LH 2 tank self pressurization testing 

• Application to Lunar Landing Vehicle (currently named Altair) 

- Parametric analysis on pressurant conditions 

- Unequal tank pressurization and draining for multiple tank designs 
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Cryogenic Propellant and Pressurization Modeling 


♦ CPPPO Model Development 

• Control Volume Definitions 

- Ullage gas: mixture of propellant 
vapor and Helium pressurant 

- Surface: infinitely thin layer 
containing no mass at propellant 
saturation conditions (used to 
pass heat between ullage and 
liquid and allow for surface layer 
evaporation) 

- Liquid: liquid propellant remaining 
in tank 

• 5 Node Model 

- Tank wall exposed to ullage gas 

- Ullage gas 

- Surface layer 

- Liquid 

- Tank wall exposed to liquid 
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Cryogenic Propellant and Pressurization Modeling 


♦ CPPPO Model Development 

• Variable Definitions 


Variable 

Units 

Definition 

W-autog 

lbm/s 

Autogenous GH2 or G02 pressurization mass flow rate 

™He 

Ibm/s 

Helium pressurization mass flow rate 

rilvent 

lbm/s 

Vent relief mass flow rate 

Qexwg 

Btu/s 

Heattransfer rate from external environment to metallictank wall skin exposed to ullage gas node 

Q«g 

Btu/s 

Heat transfer rate between ullage gas node and tank wall exposed to ullage gas node 

Q B S 

Btu/s 

Heattransfer rate between ullage gas node and saturated surface layer node 

ITlevap 

Ibm/s 

Propellant evaporation mass flow rate 

rh bb 

Ibm/s 

Propellant bulk boiling mass flow rate 

K 

Ibf-ft/s 

Ullage gas node work rate 

Qexwl 

Btu/s 

Heattransfer rate from external environment to metallictank wall skin exposed to liquid node 

Qwl 

Btu/s 

Heattransfer rate between liquid node and tank wall exposed to liquid node 

Qis 

Btu/s 

Heattransfer rate between liquid node and saturated surface layer node 

d^cond 

Ibm/s 

Propellant vapor condensation mass flow rate 

w l 

Ibf-ft/s 

Liquid node work rate 

rhi 

Ibm/s 

Propellant liquid mass flow rate 

Til TVS 

Ibm/s 

Propellant Thermodynamic Vent System (TVS) mass flow rate 

Ptg 

psia 

Tank Pressure 

% 

fl 3 

Ullage gas volume 

V, 

ft 3 

Liquid volume 


m autog 0r ™H S 
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Cryogenic Propellant and Pressurization Modeling 


♦ CPPPO Model Development 

• Various Model Information 

- Finite Difference method 

- All fluid properties updated at 
each time step using NIST 
Refprop 

- All internal tank heat transfer is 
Free Convection 

- Helium pressurant tank assumed 
isentropic blowdown 

- Tank wall is modeled as a lump 
thermal mass 

- External heat transfer rates 
assumed averaged over exposed 
tank surface area 

- Simplistic Thermodynamic Vent 
System (TVS) modeled during 
any in-space coast phases 
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Cryogenic Propellant and Pressurization Modeling 


♦ CPPPO Model Development 

• Conservation Equations 
- Ullage Gas: 


dV t d ( \ 

e, - e, * - ■ a - * k“. ) 

- Applying Ideal Gas assumption to internal energy in ullage: 


dT 

m <£v„ -^- = Q g s~ Q, g + m au,o g 


{h»,o g ~ C T tg )+ m He {h He - C T tg )+ (A eV(? , - C T tg ) 


+ Kb (*W, - T , g )- Koni { h co„d ~ C ,, g T , g )- Ken, (\ " C v, s ^ )“ 


dV w dC. 

P_ — - - m.T 


v tg 


tg 


dt 


tg tg 


dt 


- Pressure and Volume rate of change: 

p, g = pM 

K = v, +—{K b +£a P +Kvs-Kond)-v, & 

Pi Pi 
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Cryogenic Propellant and Pressurization Modeling 



♦ CPPPO Model Development 

• Conservation Equations 

- Liquid: 

dv d / x 

Qwl ~ Qh + AonAond ~ ™bb A K ~ ™bb\b ~ AvsKs ~ ~ P ,g \ m , U , ) 

at at 

- Apply internal energy and enthalpy relationship to get: 


m, 


dh } 

dt 


= QM-QlAAonA h cond- h t)-™bbA-™bb{ h bb -kVA 'Vs( h TVS - h l)+ V l 


dP. 


tg 


dt 


- Therefore, liquid propellant remains a real fluid. 

- This is very important for LH 2 during in-space coast phases since it is highly 
compressible, contracting and swelling from pressure and heat. 
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Cryogenic Propellant and Pressurization Modeling 


♦ CPPPO Model Development 

• Heat T ransfer 

- Internal Tank Heat Transfer is all in form of Free Convection 

• During rapid pressurization this assumption breaks down due to forced convective 
motion of pressurant gas in ullage 

• However, these time scales are minimal compared to the burn duration and coast 
times of interest 

- Free Convection coefficient and exponent taken from standard text with 
acceleration vector parallel and normal to plane 

• Mass T ransfer 

- Included in form of evaporation, bulk boiling, and condensation 

- Certain conditions must be achieved before any form of mass transfer occurs 

• Heat and Mass T ransfer are most difficult to validate 

- Limited data exists for pressure rise in closed tanks however the method of heat 
transfer (free/forced convection, conduction) and mass transfer (evaporation, 
bulk boiling, condensation) causing such a pressure rise is difficult to pinpoint 
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Percent Error [%] 


Cryogenic Propellant and Pressurization Modeling 


♦ CPPPO Model Validation 

• Conservation of mass and energy internal to model 

- Validated during all scenarios for heat and mass transfer 



Perr_energy % 
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Cryogenic Propellant and Pressurization Modeling 



♦ CPPPO Model Validation 

• Comparison to existing analytic model 

- ROCETS used for design of Ares I Upper Stage 

- LH 2 Tank Examples 



Time [sec] 
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Cryogenic Propellant and Pressurization Modeling 



♦ CPPPO Model Validation 


• Comparison to existing analytic model 

- ROCETS used for design of Ares I Upper Stage 

- L0 2 Tank Examples 



Cryogenic Propellant and Pressurization Modeling 


♦ CPPPO Model Validation 

• Flight Experiment 

- AS-203 (S-IVB stage orbital test, LH 2 tank results) 



-* Ptg psia -a- Pvapor psia 
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Cryogenic Propellant and Pressurization Modeling 


♦ CPPPO Model Validation 

• Ground Experiment 

- NASA Glenn LH 2 tank self pressurization testing (29% fill level) 
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Cryogenic Propellant and Pressurization Modeling 


♦ CPPPO Model Implementation 

• Applied to Altair Descent Module 

- Engine Operation Phase 


- Ptg psia 

- Pvapor psia 


-K-Ttg R 
-a-TI R 
-a-Ts R 


r ravlsdot [bmfsec 
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Temperature [R] 


Cryogenic Propellant and Pressurization Modeling 


♦ CPPPO Model Implementation 

• Applied to Altair Descent Module 
- In-Space Coast Phase 


-*^TI R 
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Cryogenic Propellant and Pressurization Modeling 



♦ CPPPO Model Implementation 

• Applied to Altair Descent Module 

- Differential Tank Draining 

- Variations in pressurization and feed line resistances will cause variations in tank 
pressures when using multiple propellant tanks 

- Result is differential draining, possibly causing early engine cutoff and CG shifts 

-*-mldot_t1 Ib/s 
-A-mldot_t2 Ib/s 
-o-mldot_t3 Ib/s 
-?-mldot_t4 Ib/s 
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Cryogenic Propellant and Pressurization Modeling 


♦ CPPPO Model Implementation 

• Other Altair Descent Module Applications 

- Used in feed system pressure loss calculations 

• Added ability to calculate pressure loss in varying lines and components 

• Also part of differential draining analysis 

- Propellant Scavenging post landing 

• Calculated various methods of venting and heating the tanks after landing to supply 
residual propellant to life support and power subsystems 

- Cold Helium tank characteristics 

• Calculated real fluid property results of rapid blowdowns of cold Helium tanks 

• Results show drastic variation from ideal fluid behavior 

- Applied to alternative tank designs 

• Used code to model pressurization characteristics of alternative and novel tank 
designs in attempt to reduce the number of propellant tanks (and avoid differential 
draining issues) 

- Liquid Methane based vehicles 

• Used model to calculate some preliminary results of pressurization for liquid methane 

- Storable vehicles 

• Have updated model to also calculate pressurization properties of storable propellant 
vehicles such as NTO/MMH 

• NIST Refprop does not contain properties for these propellants so manually coded 
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Cryogenic Propellant and Pressurization Modeling 


♦ Future Work 

• Validation against some preliminary CFD data is underway 
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♦ Questions? 
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